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We show that the daily average air humidity fluctuations exhibit non-trivial 1 //" behaviour which 
different from the spectral properties of other meteorological quantities. This feature and the 
fractal spatial structure found in clouds make it plausible to regard air humidity fluctuations as a 
manifestation of self-organized criticality. We give arguments why the dynamics in air humidity 
can be similar to those in sandpile models of SOC. 



There is widespread interest in understanding the dy- 
namic processes of large, spatially extended systems. 
There are many geophysical processes in nature where 
scale-invariance is observed. Recently, the Gutenberg- 
Richter law of earthquakes [Q, the volcanic activity ||] 
and various astrophysical phenomena have been re- 
ported to show scale-invariant properties. However, such 
quantitative analysis is often limited by a shortage of 
available data or a lack of understanding of the basic 
effects. One of the largest extended systems is our atmo- 
sphere. Although all the basic mechanisms that govern 
the dynamics of the atmosphere have been well known for 
quite a long time, a detailed understanding and adequate 
characterization of the fluctuations of various statistical 
quantities of the lower atmospheric boundary layer are 
still not complete. The conventional viewpoint is that the 
underlying mechanism is random and can be considered 
as an autoregressive process. Quite surprisingly, in Ref. 
Q it was possible to demonstrate that the distribution 
and spectral properties of the daily medium temperature 
fluctuations coincide with those measured in a helium cell 
in Refs. |^,^. The daily mean temperature fluctuations 
on a one day scale and the Soft Turbulent 2K helium cell 
fluctuations on a millisecond scale look the same since 
the dimensionless parameters, the Prandtl and Raylegh 
numbers of the two systems, are the same. The main 
objection to using daily average data for deeper anal- 
ysis is the belief that random weather changes (fronts, 
etc.) dominate the temporal behaviour of the quantities 
measured at a typical meteorological station. In real- 
ity, in the case of daily temperature fluctuations, it was 
observed that the meteorological fronts due to the large 
scale coherent spatial motion affect only the short (1-10 
days) timescale, while the vast majority of timescales (10- 
10000 days) is governed by the inherent fluctuations of 
the Benard- Raylegh convection in the lower 500m bound- 
ary layer H]. In this Letter we would like to show that 



the fluctuations of the daily average relative air humid- 
ity data also reflect some inherent physical processes and 
that this measurement is an ideal candidate for observing 
large scale-invariance in the atmosphere. 

The air humidity can be quantified by the partial va- 
por pressure py of gaseous water in the air. At a given 
temperature T, the amount of water in the air is limited. 
The maximum partial vapor pressure corresponding to 
the maximum water content is the saturation vapor pres- 
sure Psat(T). The relative air humidity e at temperature 
T is the ratio of the actual to the saturation pressure 
measured in percent 

e = -^jTfT, ■ 100%. (1) 

Psat[T) 

A detailed summary of the measurement and the defini- 
tion of the air humidity and the saturation vapor pressure 
are given in the classic work of Penman Q . 

The special importance of the humidity from the point 
of view of self-organized criticality (SOC) § lies in its 
threshold behaviour. If the actual humidity in the air lo- 
cally exceeds 100%, the water starts condensing. In the 
absence of macroscopic spatial motion of the air, the mi- 
croscopic amount of condensed water is transported on 
microscopic scales by diffusion. It can reach a new lo- 
cation, where the local relative humidity is under 100% 
and can return to the gaseous state. If the humidity 
in the neighbourhood is also saturated, the amount of 
condensed water will increase. This can initiate a chain 
reaction like the avalanches in the 'sandpile' models 
We propose that large decreases in the air humidity are 
consequences of such processes. On the other hand, the 
humidity in the air is coming from open water surfaces, 
from the moisture in soil and from plants and animals. 
If the air humidity is under 100%, the evaporation from 
these sources continually increases the water content of 
the air. This is the analog of dropping sand on the sand- 
pile, and this drives the system toward the fully saturated 
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100% humidity state which is the stationary and critical 
state of the system at the same time. It is very likely 
that these aspects of the dynamics help to build up a 
self-organized critical state on large scales. The two main 
characteristics of the SOC behaviour is the 1//" type be- 
haviour of the power spectra of the fluctuations, and the 
spatial scaling behaviour manifested in the formation of 
fractal structures |8|. One of the most well known frac- 
tals in nature are clouds which are formed in the 
chain reaction-like condensation of air humidity. These 
processes and large scale fractal behaviour of clouds can 
be observed in meteorological satellite pictures |jl^ . The 
condensation of air humidity into clouds can be modeled 
by cellular automata, which has recently been reported 
to show SOC behaviour [pT| . 

To detect scaling in the temporal behaviour, we have 
analyzed relative air humidity data collected by the Hun- 
garian Meteorological Service in various stations in Hun- 
gary during the period 1963 - 1988. The main difhculty in 
analyzing meteorological time series is that seasonal peri- 
odic variations can cause peaks in the Fourier transforms 
and make it difficult to determine the power spectra of 
the fluctuations. Figure 1 shows the seasonal variation 
e{i) of a typical station (Balaton, county Heves), which 
is defined as 

^ 1988 

^(*)=26 ^ ""^'^y^' i = l,...,365 (2) 

a=1963 

where e{i,y) denotes the daily average air humidity of the 
ith day of the yth year measured at the meteorological 
station. To diminish the effect of seasonal variations we 
subtract it from the daily data. Then the fluctuation e/ 
is given by the deviation of the actual humidity from the 
seasonal variation ^ : 

ef{i,y) = e{i,y) - e{i), (3) 

where the notation is the same as above. In this way we 
obtain a 26 x 365 day long record of daily average relative 
air humidity fluctuations. Hungary can be considered as 
meteorologically homogeneous from the point of view of 
fluctuations |Q due to the geography. Different stations 
measure coherent fluctuations Q, and one station can 
be considered as representative. Figure 2 shows a typi- 
cal record of the daily average humidity fluctuations for 
the year 1969. Figure 3 shows a histogram of the daily 
average humidity fluctuations, which is clearly a Gaus- 
sian distribution, with zero mean and standard devia- 
tion of CT = 9.34%. The power spectrum 5'(/), which is 
the squared magnitude of the Fourier amplitude, can be 
obtained by standard numerical methods |l^. Figure 4 
shows the power spectrum of the time series measured in 
Balaton. Almost the whole frequency range can be fitted 
by a simple power law 

Sif) - ^, (4) 



with a = 0.61 ± 0.01. In addition to the large scale spa- 
tial scaling behaviour this shows the existence of scaling 
behaviour also in the time domain. A very nice feature of 
this result is that meteorological fronts and other uncor- 
related events do not seriously affect the spectra. If the 
underlying process is SOC, the exponent a can be related 
to the scaling of the weighted average of the avalanches 
[Q . The weighted average is defined as 

A(t)=^,s^P(.,t), (5) 

s 

where P(s, t) is the distribution of the avalanche sizes s 
as a function of the duration t of the avalanches. For 
SOC this scales as 

m ~ ^^ (6) 

where the exponent satisfies — 1 < fi < +1. The 
exponents fi and a are related |14| as a ^ fi + I, which 
yields n = -0.39 ±0.01. 

Finally, to demonstrate the non-trivial nature of a 
1//" type power-density spectrum, on Figures 5 and 6 
we show the power spectra of the daily temperature and 
rain fluctuations respectively, measured at the same me- 
teorological station. The temperature record has been 
analyzed in Rcf. while the rain record appears to be 
white noise (ie. a constant power spectrum). Neither 
show any sign of scaling behaviour. 

In summary, we have found that the spectral density 
of daily average air humidity fluctuations shows scaling 
behaviour. This is very different from the spectral prop- 
erties of other meteorological quantities. This and the 
fractal behaviour of the condensed air humidity in space 
make it very plausible to regard air humidity fluctuations 
as a manifestation of self-organized criticality. We gave 
arguments why the dynamics in air humidity can be simi- 
lar to those in sandpile models of SOC. We hope that this 
study can inspire more detailed analysis of the space-time 
dynamics of clouds and air humidity fluctuations. 
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FIG. 1. The seasonal variation of the daily average air hu- 
midity (Eq.(2)) obtained from the data of the meteorologi- 
cal station Balaton, county Heves, Hungary for the period 
1963-1988. 



FIG. 2. Air humidity fluctuations (Eq.(3)) measured in 
1969. 



FIG. 3. Unnormalized histogram of the air humidity fluc- 
tuations measured in the period 1963-1988. The dashed line 
is the fitted Gaussian distribution with standard deviation 
a = 9.34% 

FIG. 4. Unnormalized power spectrum of the daily average 
relative air humidity fluctuations. The dashed line is the best 
fit given by (Eq.(4)). 



FIG. 5. Unnormalized power spectrum of the daily medium 
temperature fiuctuations. 

FIG. 6. Unnormalized power-density spectrum of the daily 
rain fluctuations. 
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